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In this paper, by making use of Duan's topological current theory, the evolution of the vortex 
filaments in excitable media is discussed in detail. The vortex filaments are found generating or 
annihilating at the limit points and encountering, splitting, or merging at the bifurcation points of 
a complex function Z{x, t). It is also shown that the Hopf invariant of knotted scroll wave filaments 
is preserved in the branch processes (splitting, merging, or encountering) during the evolution of 
these knotted scroll wave filaments. Furthermore, it also revealed that the "exclusion principle" in 
some chemical media is just the special case of the Hopf invariant constraint, and during the branch 
processes the "exclusion principle" is also protected by topology. 

PACS numbers: 02. 10. Kn, 82. 40. Ck, 02. 40. Xx, 03. 65. Vf 



I. INTRODUCTION AND MOTIVATION 

Scroll waves are three-dimensional (3D) extensions of 
the familiar spiral waves of excitable media. They have 
been observed in a variety of physical, chemical and bi- 
ological systems[l|, 0, S S S @- Recently, scroll wave 
have drawn great interest due to its importance in the 
the mechanism of some re-entrant cardiac arrhythmias 
and fibrillation which is the leading cause of death in 
the industrialized world [3, 0, @]- The scroll wave rotate 
about a linelike filaments called vortex filament, and usu- 
ally can be defined in terms of a phase singularity. In 
three-dimensional excitable media, the vortex filament is 
commonly a closed ring, and these vortex filaments can 
form linked and knotted rin gs w hich contract to compact, 
particle-like bundles 0, [HI, El, [H [ll E H [13, E E 

Control of scroll wave is a more important and complex 
problem for all excitable media. The dynamics of a 3D 
scroll wave are determined not only by properties of the 
excitable media but also by the geometry and topology 
of vortex filament [g, 0, H, ISj]. This implies that the con- 
trol of scroll wave should stronger depend on our under- 
standing on scroll wave topology. This inspirits us to use 
the topological viewpoint to s tudy the scroll wave topol- 
ogy. In previous works [ll,[li,[l3,[ll[li,[2SL21, 22, 2j], 
many authors have made great contributions to this issue 
and employed topological arguments to understand scroll 
wave. Some most important topological constraints on 
behaviors of the vortex filaments have been investigated. 
These topological rules may have some important appli- 
cations in practice. In particular, the topological con- 
straint on knotted vortex filaments is believed to relate 
to topological characteristic numbers of knotted vortex 
filament family, such as the winding, the self-linking and 



the linking numbers. In Ref.[19|], Winfree and Strogatz 
have proposed an "exclusion principle" which governed 
the scroll wave knotting and linking through each other 
in chemical system. This exclusion principle gives the 
constraint on the linking numbers and winding numbers 
of scroll wave. 
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Recently, Duan's topological current theorv[15lll6l.l23l. 
Y2A . l25l . l2a | has been applied to study the topological 
properties of spiral waves and scroll waves. Zhang et 
al. ilal presented a rigorous topological description of spi- 
ral waves and scroll waves. They derived precise expres- 
sions of spiral wave and scroll wave topological charge 
density. Based on their work, we study the branch pro- 
cess of the spiral waves and calculate the knotted in- 
variant for knotted vortex filaments by using the Duan's 
topological current theory, we proposed that the knot- 
ted invariant (which is just the Hopf invariant) may im- 
ply a new topological constraint on scroll wave lm [23| . 
However, in Ref.[l5| and in our previous work[l6[, the 
discussions are based on an important regular condition 
D{-) ^ 0. When this condition fails, what will happen? 
Main purpose of this paper is to detail this problem. 



In this paper, by making using Duan's topological the- 
ory, firstly, we will extend our branch theory of spiral 
wave in 2D to 3D scroll wave, and study the generating, 
annihilating, colliding, splitting and merging of vortex 
filaments from a topology viewpoint. Secondly, based on 
the branch process of vortex filaments, it is showed that 
the Hopf invariant of knotted scroll wave is preserved in 
the branch process. This is consistent with our proposal 
that the Hopf invariant may implies a new constraint on 
scroll wave. Third, it is also shown that the "exclusion 
principle" in some chemical media is just the special case 
of the Hopf invariant constraint, and during the branch 
processes the "exclusion principle" is also protected by 
topology. 



II. TOPOLOGICAL STRUCTURES OF VORTEX 
FILAMENTS 

In order to maintain the continuity of the whole work 
and make the background of this paper clear, in this sec- 
tion, we give a brief review of the topological current the- 
ory of vortex filaments. We chose to work with a general 
two- variable reaction-diffusion system whose mathemati- 
cal description in terms of a nonlinear partial differential 
equation. This equation is written as 



dtu = /(u,w) -I- £)„V^M, 
dtv = g{u,v) + Dy\/'^v, 



(1) 



where u and i; represent the concentrations of the 
reagents; V'^ is the Laplacian operator in three- 
dimensional space; f{u,v) and g{u,v) are the reaction 
functions. Following the description in Ref . [lai. [27| . we 
define a complex function Z — (ji^+icj)'^, where (j)^ = u—u* 
and 0^ = V — V* . Here u* and v* are the concentrations 
of the vortex filaments. 

We know that the complex function Z = 0^ + icj)^ 
can be regarded as the complex representation of a two- 
dimensional vector field Z = {(}?-, (fp'). Let us define the 
unit vector: n'' = |^(a = 1,2; ||0f = (t)°'(t)°' = Z*Z). It 
is easy to see that the zeros of Z are just the singularities 
of n. Using this unit vector n, an "induced abelian gauge 
potential" can be constructed with 



A^, = eabn'^d^n", ^^ = 0,1,2,3; 



(9o,V), do^dt 



dt' 



(2) 



the gauge field strength given by this gauge potential is 



-t' ^1/ — o^Au — OijA^. 



(3) 



According to Ref. [13, M, S, H, 123, 111 , the two dimen- 
sional topological tensor current is defined as 



1 



i.-^-.^ -.V 



,M^Ap^ 



1 

2^" 



.M^Apn „a^ &^ 



(4) 



It is easy to show that the topological tensor current K^'^ 
can be rewritten in a compact form 



X 

where D^'^ is the general Jacobian determinants 



(5) 



where -D'(f ) = ^^^ £abdj4>°'dk(f> is the Jacobian vector. 
This delta function expression of the topological current 
j' tells us it doesn't vanish only when the vortex filaments 
exist, i.e.. 



] 



= 0, if and only if 7^ 0; 
^ 0, if and only ii (j> — 0. 



(9) 



So the sites of the vortex filaments determine the nonzero 
solutions of j^ . The implicit function theory shows that 
under the regular condition|28| 

D^-{^)^Q, (10) 

X 

the general solutions of 

0^(t,x)=O, (j)^{t,x)^0 (11) 

can be expressed as 

x^^xUt,s), x^^xf{t,s), x^=xf{t,s), (12) 

which represent the world surface of N moving isolated 
vortex filaments with string parameter s (l = 1,2, •,7V). 
These singular strings solutions are just the vortex fila- 
ments. In delta function theory [29|, one can prove that 
in three-dimensional space. 






1=1 



(13) 



where D(^) = ieJ^,„„M^Ml and S, is the l-th planar 
element transverse to Li with local coordinates (u^,u^). 
The positive integer Pi is the Hopf index of (/)-mapping, 
which means that when x covers the neighborhood of 
the zero point xi{s,t) once, the vector field (f> covers the 
corresponding region in (j) space for Pi times. Meanwhile 
the direction vector of Li is given by[23, 123] 



dx' 
ds 






(14) 



Then considering Egs. fTS)) and Eqs.jM]), we obtain the 
inner structure of/'. 



f = <5(0)i?H-) 



N 

= ^Pim / dx'-d^ix-xi), 



(15) 



e''''Df"'{-) ^ e'^'^Pdxcl)'' 



Defining the spatial components of K'^" as 



(6) 



1 

2^' 



/ = K°' ^ —^Jk^^^Q^^-d^n\ t,j, k = l,2, 3, (7) 



we have 



x 



(8) 



where 77; — sgnD{-) = ±1 is the Brouwer degree of 
0-mapping, with r]i — I corresponding to the vortex fil- 
ament and rji — —I corresponding to the antivortex fila- 
ment. We find that the topological current j is just the 
charge density vector p of the vortex filament in Ref. [l3] ■ 
In our theory, the topological charge of the vortex fila- 
ment Li is 



U = 



= j-da^Wi ^ Pim, 



(16) 



in which Wi is just the winding number of (j> around L;, 
the above expression reveals distinctly that the topo- 
logical charge of vortex filament is not only the wind- 
ing number, but also expressed by the Hopf indices and 
Brouwer degrees. The topological inner structure showed 
in Eq. ()16|) is more essential than usually considered and it 
will be helpful as a complement of the current description 
of scroll wave only by winding number in topology. This 
shows that the advantage of our topological description 
of the vortex filaments. 

In the above, we give an prime introduction of the 
topological current theory of scroll wave, which is based 
on the rigorous mathematics background that we called 
Duan's topological current theory. Here we give some re- 
marks about above results in order, (i). The definition of 
the complex function Z = (j)^ + icj)^ originates from u and 
V. Therefore, Z is the function that shows the concentra- 
tion distribution of the reagents. The exact expression 
of Z is determined by the nonlinear partial differential 
equation ^ . If the exact analytical solution of Eq. ([T|) is 
known, we can directly calculate the zero points of Z and 
determine the Hopf indices and Brouwer degrees, which 
describe the inner topological structures of the vortex fil- 
aments of scroll wave. (ii). The traditional theoretical 
description of scroll wave usually used the phase singular- 
ity method and cannot directly deal with the zero points 
of Z . The above results show that the topological cur- 
rent theory of scroll wave provide a available approach 
to deal with the zero points, (in). The regular condition 
(fTO|) will plays a essential role in determining the stabil- 
ity of the scroll wave. This condition also depends on 
the exact solution of the Eq.([T]). It can be determined 
when the expression of Z is known. How this condition 
determine the stability of the scroll wave is the tasks of 
the following sections. 



III. THE BRANCH PROCESS OF VORTEX 
FILAMENTS AT THE LIMITED POINT 

However, from the above discussion we know that 
the results mentioned are obtained under the condition 
£)/ji'(£) _^ Q^ When this condition fails, i.e., the Brouwer 
degrees rji are indefinite, what will happen? In what fol- 
lows, we will study the case when D'^'^{-) = 0. It often 

happens when the zero of Z includes some branch points, 
which lead to the bifurcation of the topological current. 
Generally speaking, the evolution of a vortex filament 
Li can be discussed from Eq.®. From Eq.(|4]), consid- 
ering that e^'^^'' is a fully antisymmetric tensor, we can 
prove that 



d^Kf''' = 0, 



that is 



do.f 



djK^' 



0. 



(17) 



(18) 



This is the continuity equation constraints on vortex fil- 
aments. In order to discuss the evolution of these vortex 



filaments and simplify our study, we fixed the x^ = z 
coordinate and take the XOY plane as the cross section, 
so the intersection line between the Li 's evolution surface 
and the cross section is just the motion curve of L/. In 
this case the 2D topological current is defined as 



f = X"3 ^ s^(^cb)D'>{^) 



and 



K' = K'- 



5\c^)D\^, i = l,2. 



(19) 



(20) 



It is easy to see that j^ and K^ satisfy the continuity 
equation 



dof 



d,K' = 0. 



(21) 



The velocity of the intersection point of vortex filament 
and the cross section is given by 



dx' 



dH4>/x) 



D°{(j)/x)' 

From Eq. (f22|) it is obvious that when 
D^{(j)/x) = 
at the very point (i*, x*), the velocity 

dx'^ _ D'^{(l)/x) 



(22) 



dx^ _ D^{4)/x) 
~dt ^ 



D^{<t)/2 



D^{4)/x) 



(23) 



(t*,x*) 



is not uniquely determined in the neighborhood of 
{t*,x*). This critical point is called the branch point. 
In Duan's topological current theory usually there are 
two kinds of branch points, namely the limit points and 
the bifurcation points, each kind of which corresponds to 
different cases of branch process. 

First, in this section, we only study the case that 
the zeros of the complex function Z includes some limit 
points which satisfy 



D\^) 



^0, DH^) 



(t*,s-) 



^0. 



(24) 



(t*,s*) 



For simpficity, we assume that D'^{(j)/x)\(^f_s*) ^ 
is always satisfied in our next discussions. When 
D^{(t)/x)\(^t-',x-') ^ 0, from Eq.flSl) we obtain 



dx^ 

~dr 



D\<f>/x) 



DO{(l)/x) 



{t',x') 



dt 

d^ 



= 0. 



(25) 



(26) 



{f.x') 



Then, the Taylor expansion of of t = t{x^ , t) at the limit 
point (i*,x*) of vortex filament, one can obtain 



t-t* 



1 dH 



2 (da;i) 



(xi 



(27) 



(t*,s*) 



which is a parabola in x^ ~ t plane. From Eq. (|27|) we 
can obtain two solutions x\{t) and X2{t), which give two 
branch solutions of vortex filament at the limit points. If 



dH 



(dx^) 



1\2 



>o, 



(28) 



(t%x*) 



we have the branch solutions for t > t* [see Fig. 1(a)]; 
otherwise, we have the branch solutions for t < t*[see 
Fig. 1(b)]. The former is related to the origin of the 
vortex filament at the limit points, and the later is the 
annihilation of the vortex filament. At the neighborhood 
of the limit point, we denote the length scale I = Ax. 
From Eq.(|27p, one can obtain the approximation relation 



1<X\\ t~t* ||l/2 



(29) 



The growth rate 7 = 1/ Ax or annihilation rate of vortex 
lines is 



7(x(i-r)-i/^ 



(30) 



From the continuity equation Eq. (|2ip . we know that 
the topological number of the vortex filament is identi- 
cally conserved. This means that the total topological 
number of the final vortex filaments equals that of the 
initial vortex filaments. The total numbers of these two 
generated or annihilated vortex filaments must be zero 
at the limit point; i.e., the two generated or annihilated 
vortex filaments have be opposite. 



Pim + (32f]2 = 0, 



(31) 



which shows that (3i = (32 and rji = — ?72- One can see 
the fact that the Brouwer degree 77 is indefinite at the 
limit points implies that it can change discontinuously at 
limit points. 

For a limit point it is required that £'^(0/2:)|(f,s*) 7^ 
0. As to a bifurcation point [33|, it must satisfy a more 
complex condition. This case will be discussed in the 
following section. 



IV. THE BRANCH PROCESS OF VORTEX 
FILAMENTS AT THE BIFURCATION POINT 

Now let us study the bifurcation of the vortex line at 
its bifurcation point where 



^"0 



0, D\^) 



{t',x'} 



0. 



(32) 



(t',x') 



These two restrictive conditions will lead to an important 
fact that the function relationship between t and x^ is 
not unique in the neighborhood of the bifurcation point 
{t*,x*). The equation 



(hP__ D\(f>/x) 
~dt ~ mjjx) 



(33) 



(a) 



FIG. 1: We fixed the limit point {x^*,t*) at the origin of 
{x^ — t) plane, (a) The branch solutions for Ea. (|27p when 
d?t/ {dx^Y\(^t* ^s*) > 0, i.e., a pair of vortex filaments with 
opposite charges generate at the limit point, i.e., the origin of 
vortex filaments, (b) The branch solutions for Eq. (|27|) when 
d^t/ {dx^)'^\(^t* ,x*) < 0, i.e., a pair of vortex filaments with 
opposite charges annihilate at the limit point. 



which, under restraint of Eg. ([5^ . directly shows that the 
direction of the integral curve of Eq. ([33]) is indefinite at 
the point {t*,x*). This is why the very point (t*,a?*) is 
called a bifurcation point. 

Assume that the bifurcation point {t* , x* ) has been 
found from Eqs. (fTT|) and ((32|) . We know that, at the 
bifurcation point (t*,a;*), the rank of the Jacobian ma- 
trix [d(t>/dx\ is 1. In addition, according to the Duan's 
topological current theory, the Taylor expansion of the 
solution of Ea. fTTjl in the neighborhood of the bifurca- 
tion point {t*,x*) can be expressed as|23l[25| 

A{x^ - x*^f + 2B{x' - x*^){t - r ) + C{t - t*f = Qf34) 



which leads to 



dT^ 



dx^ 



and 



^(:^)^ + 2i.(:^: 



'dx^ 



'dx^ 



A = 0, 



(35) 



(36) 



(t*,s*) 



where A, B, and C are three constants. The solutions of 
Eq. ([55)1 or Eq. ([35]) give different directions of the branch 



curves at the bifurcation point. There are four possi- 
ble cases, which will show the physical meanings of the 
bifurcation points. 




FIG. 2: We fixed the bifurcation point (a;^*,i*) at the origin 
of (x^ — t) plane. Two vortex filaments meet and then depart 
at the bifurcation point. 

Case 1 (A ^ Q). For A = A{B^ - AC) > from 
Eg. ([55)) we get two different directions of the velocity 
field of vortex filaments 



dt 



B ± VB^ - AC 



(37) 



(t*,s*) 



Which is shown in Fig. 2. It is the intersection of two vor- 
tex filaments with different directions at the bifurcation 
point, which means that two vortex filaments meet and 
then depart from each other at the bifurcation point. 

Case 2 (A ^ 0). For A = 4(5^ - AC) = from 
Eg. ([35]) we obtain only one direction of the velocity of 
vortex filaments 



djp_ 
dt 



(**,£*) 



-B 



(38) 



which includes three important situations. (a)Two vor- 
tex filaments tangentially encounter at the bifurcation 
point [See Fig. 3(a)]. (b)Two vortex filaments merge 
into one vortex filament at the bifurcation point [See 
Fig. 3(b)]. (c) One vortex filament splits into two vor- 
tex filaments at the bifurcation point [See Fig. 3(c)]. 

Case 3 (A = 0,C ^ 0). For A = 4(5^ - AC) = 0, we 
have 



dt 
dx^ 



-B±VB^-AC 
C 



0, 



2B 

c ■ 



(39) 



There are two important cases: (a) One vortex filament 
splits into three vortex filaments at the bifurcation point 
[See Fig. 4(a)]. (b) Three vortex filaments merge into one 
vortex filament at the bifurcation point [See Fig. 4(b)]. 

Case 4 (A = C = 0). EquationdSS]) and Eqjlll) give 
respectively 



dx^ 
~dF 



0, 



dt 
d^ 



0. 



(40) 



(a) 






FIG. 3: We fixed the bifurcation point {x *, i*) at the origin of 
(x^—t) plane, (a) Two vortex filaments tangentially encounter 
at the bifurcation point, (b) Two vortex filaments merge into 
one vortex filament at the bifurcation point. (c)One vortex 
filament splits into two vortex filaments at the bifurcation 
point. 



This case is obvious similar to Case 3, see Fig. 5. 

The above solutions reveal the evolution of the vortex 
filaments. Besides the encountering of the vortex fila- 
ments, i.e., a vortex filament pair encounter and then 
depart at the bifurcation point along different branch 
cures [See Fig. 2 and Fig. 3(a)], it also includes splitting 
and merging of vortex filaments. When a multi-charged 
vortex filament moves through the bifurcation point, it 
may split into several vortex filaments along different 
branch curves [See Fig. 3(c), Fig. 4(a) and Fig. 5(b)]. On 
the contrary, several vortex filaments can merge into a 





FIG. 4: We fixed the bifurcation point (a;^*,i*) at the origin 
of {x^ — t) plane, (a) One vortex filament splits into three 
vortex filaments at the bifurcation point, (b) Three vortex 
filaments merge into one vortex filament at the bifurcation 
point. 



vortex filament at the bifurcation point [See Fig. 3(b) and 
Fig.4(b)]. 

At the neighborhood of the bifurcation point, we de- 
note scale length I = Ax. From Eqs. ([57| - (|5^ we can 
then obtain the approximation asymptotic relation 



;cx {t-t*). 



(41) 



The growth rate 7 or annihilation rate of vortex filament 
7 of the vortex filament is 



7 oc const. 
From Eg. (1^0)1 ■ one can obtain 

I — const, 7 = 0. 



(42) 



(43) 



It is obvious that the vortex filaments are relatively at 
rest when I — const. 

The identical conversation of the topological charge 
shows the sum of the topological charge of these final 
vortex filaments must be equal to that of the original 
vortex filaments at the bifurcation point, i.e., 



E/^'>^'> = E/3// 



Vif 



(44) 



FIG. 5: We fixed the bifurcation point (a;^*,i*) at the origin 
of {x^ — t) plane. In this case, Two vortex filaments intersect 
normally at the bifurcation point, (a) Three vortex filaments 
merge into one vortex filament at the bifurcation point, (b) 
One vortex filament splits into three vortex filaments at the 
bifurcation point. 



for fixed /. Furthermore, from the above studies, we see 
that the generation, annihilation, and bifurcation of vor- 
tex filaments are not gradually changed, but suddenly 
changed at the critical points. 



V. HOPF INVARIANT CONSTRAINT ON 
SCROLL WAVE 

In this section, we will research the topological prop- 
erties of the knotted vortex filaments. We first consider 
the continuity equation constraint, from Eq.([7]) one can 
obtain 



^^f = 0, 



(45) 



which can also be derived from Ea. (|17p . The continuity 
equation (|45p implies that the vortex filament may be 
either closed loops or infinite curves. In Ref.[l5|, zhang 
et.al. pointed that the continuity equation (j45p is con- 
sistent with the topological rule which governs the scroll 
wave pinning to an inclusion 17]. 

Except the continuity equation constraint on scroll 
wave, the complex scroll wave topology may provide 



other topological requirements on scroll wave. In the fol- 
lowing discussions in this section, we will study an impor- 
tant knotted invariant which constraints on scroll wave. 
It is well know that the Hopf invariant is an important 
topological invariant to describe the topological charac- 
teristics of the knot family. In our topological theory 
of knotted vortex filaments, the Hopf invariant relates 
to the topological characteristic numbers of the knotted 
vortex filaments family. In a closed three-manifold ill 
the Hopf invariant is defined as|16l[24| 



H = ^ f AAF=^ 
27r J]^j ZTT 



A,fd^x. 



M 



Substituting Eg. pSl) into Eq. p6)l . one can obtain 



1 ^ 
1=1 



A.dx'- 



(46) 



(47) 



It can be seen that when these N vortex filaments are N 
closed curves, i.e., a family of N knots £,i{l — 1, 2, • • •, N), 
Eq.gT]) leads to 



1 



N 



H=—Y.Wif Adx\ 
"" 1=1 •^?' 



(48) 



This is a very important expression. Consider a trans- 
formation of complex function Z = e^^Z, this gives the 
U(l) gauge transformation of Ai : A,^ — Ai + di9 , where 
9 G R is a phase factor denoting the U(l) gauge trans- 
formation. It is seen that the di9 term in Ea. (|15)) con- 
tributes nothing to the integral H when the vortex fila- 
ments are closed, hence the expression ((48]) is invariant 
under the U(l) gauge transformation. As pointed out in 
Ref. [ia| , a singular vortex filament is either closed ring 
or infinite curve, therefore we conclude that the Hopf 
invariant is a spontaneous topological invariant for the 
vortex filaments in excitable media. 

According to our previous work in Ref.[l6|, a precise 
expression of the Hopf invariant is 



H 



N 

E 

k=l 



N 



W^SLi^k) 



J2 WkWiLki^k,Ci) 
k,i=i(k^i) 



(49) 



where Lk{^k,S.i) is the Gauss linking number between 
different knotted vortex filaments ^k and ^/, and SL(^k) 
is the self-linking number of closed filament ^k with an 
imaginary closed filament infinitesimally nearby 130. 1 3X1 . 
|32] . The Eq. (^^ reveals the relationship between H and 
the self-linking and the linking numbers of the vortex fila- 
ments knots family. Since the self-linking and the linking 
numbers are both the invariant characteristic numbers of 
the vortex filaments knots family in topology, H is an 
important topological invariant required to describe the 
linked vortex filaments in excitable media. 

In the following we will discuss the conservation of the 
Hopf invariant in the branch processes of knotted fila- 
ments. In the branch process of vortex filament, we note 



that the sum of the topological charges of final vortex 
filaments must be equal to that of the initial vortex fila- 
ments at the bifurcation point. This conclusion is always 
valid because it is in topological level. So we have, 

(a) for the case that one filament ^ split into two fila- 
ments ^1 and ^2 , we have W^ = W^^ + W^^ ; 

(b) two vortex filaments ^i and ^2 merge into one fila- 
ments: W^, + W52 = W^; 

(c) two vortex filaments ^1 and ^2 meet, then depart as 
other two filaments ^3 and ^4: W^^ + W^^ — ^?3 + ^U- 

In the following we will show that when the branch 
processes of knotted vortex filaments occur as above, the 
Hopf invariant is preserved: 

(A) The splitting case. We consider one knot ^ split 
into two knots ^1 and ^2 which are of the same seif-linking 
number as ^ {SL{^) = SL{^i) = SL{^2))- And then we 
will compare the two number H^ and iJ^^+^g (where H^ 
is the contribution of ^ to iJ before splitting, and H^-^^^^ 
is the total contribution of ^1 and £,2 to H after splitting. 
First, from the above text we have W^ = W^^ + W^^ in 
the sphtting process. Second, on the one hand, noticing 
that in the neighborhood of bifurcation point, ^1 and ^2 
are infinitesimally displace from each other; on the other 
hand, for a knot ^ its self-linking number SL{£) is defined 



SL{O^Lk{£,£v), 



(50) 



where £v is another knot obtained by infinitesimally dis- 
placing ^ in the normal direction y[3J]. Therefore 

SL{0 = 5L(a) = SL{£2) = ifc(ei, 6) = ifc(6, ei)(51) 
and 

Lkitek)^Lki£i,ek)^Lk{^2,ek) m 

(where ^^ denotes another arbitrary knot in the 
family(^^ 7^ ^,^^ ^ ^1.2))- Then, third, we can compare 
H^ and H^-^+^^ before splitting, 

N 

H^ = W^SL{0+ Y. 2W^W^,^Lk{tQ, (53) 

where Lfc(^,^^) = Lfc(^^,^); after splitting, 

^6+6 = WlSLi£i) + WlSLi£2) + 2W^,W^,Lk{£i,£2) 

N 

+ J2 2W^.W^,^Lk{£^,ek) 

k=l{e,^ii,2) 

N 

+ Y. 2VF5,Ty^,LMC2,CD- (54) 

k=l{e,^ii,2) 



Comparing Eqs. ([55|) and ([5i)) . we have 



■^e - -^ei+?2 



(55) 



This means that in the splitting process the Hopf invari- 
ant is conserved. 



8 



(B) The mergence case. We consider two knots ^i and 
^2, which are of the same self-linking number, merge into 
one knot ^ which is of the same self-linking number as ^i 
and ^2 • This is obviously the inverse process of the above 
splitting case, therefore we have 



H> 



C1+C2 



H, 



(56) 



(C) The intersection case. This case is related to the 
collision of two knots, we consider that two knots ^1 and 
^2, which are of the same self-linking number, meet, and 
then depart as other two knots ^3 and ^4 which are of the 
same self-linking number as ^1 and ^2- This process can 
be identified to two sub-processes: ^1 and ^2 merge into 
one knot ^, and then ^ split into ^3 and ^4. Therefore, 
from the above two cases (B) and (A) we have 



-f^Cl+C2 - ^«3+C4 



(57) 



Therefore we acquire the result that, in the branch pro- 
cesses during the evolution of knotted vortex filaments 
(splitting, mergence, and intersection), the Hopf invari- 
ant is preserved. 

The above analysis show that the branch processes of 
knotted vortex filament family must satisfy the Hopf in- 
variant constraint. This conclusion is obtained only from 
the viewpoint of topology without using any particular 
models or hypothesis. Therefore, the Hopf invariant is a 
more extensive topological constraint on scroll wave, and 
it is valid in almost systems which support the existence 
of scroll wave. 

According to Winfree and Stogatz[l9|. there is an "ex- 
clusion principle" governed the scroll wave knotting and 
linking through each other in chemical system. The 
chemical requirement plays a crucial role in such "ex- 
clusion principle". It states that the topology of scroll 
wave in such system must satisfy the constraints that: 



N 



WkSLia) 



=i(/#fe) 



6) = o. 



(58) 



It is very easy to see that the "exclusion principle" makes 
the Hopf invariant trivially, i.e., H = 0. This is just 
a special case of the Hopf invariant constraint. When 
branch processes of scroll wave occur, it is obvious from 
above discussion that the "exclusion principle" is also 
protected by topology. 



VI. CONCLUSION AND DISCUSSIONS 

First, we give a prime review of the topological theory 
of vortex filaments in three dimensional excitable media. 
When D{(j)/x) — 0, the intersection, splitting, and merg- 
ing of line defects in three-dimensional space are inves- 
tigated in detail by making using of Duan's topological 
current theory. Second, the evolution of vortex filaments 
in (3 -|- l)-dimensional space-time is studied. There ex- 
ist crucial cases of branch processes in the evolution of 



vortex filaments when the Jacobian D{(j)/x) = 0, i.e., rji 
is indefinite. At one of the limit points of the complex 
function Z, a pair of vortex filaments with opposite topo- 
logical charge can be annihilated or generated. At one of 
the bifurcation points of Z, a vortex filament with topo- 
logical charge W^ may split into several vortex filaments 
(total topological charges is W^); conversely, several vor- 
tex filaments (total topological charges is W^) can merge 
into one vortex filament with a topological charge W^. 
Also, at one of the bifurcation points of Z, two filaments 
meet and then depart. These show that vortex filaments 
are unstable at these branch points of Z. From the topo- 
logical properties of the complex function Z, we obtained 
that the velocity of the vortex filaments is infinite when 
they are being annihilated or generated, which agrees 
with the similar results of line defects what was obtained 
by Bray and MazenkopJSj. The velocity of the vortex fil- 
ament at the limit point or bifurcation point has been 
shown clearly in Fig.l to Fig. 5 (The slope of the curve). 
Third, based on the branch process of vortex filaments, 
it is showed detailed that the Hopf invariant of knot- 
ted scroll wave is preserved in the branch process, this is 
consistent with our proposal that the Hopf invariant may 
implies a new constraint (the Hopf invariant constraint) 
on scroll wave. Furthermore, it also revealed that the 
"exclusion principle" in some chemical media is just the 
special case of the Hopf invariant constraint, and dur- 
ing the branch processes the "exclusion principle" is also 
protected by topology. Finally, we would like to point 
out that both branch theory of vortex filaments and the 
Hopf invariant constraint in this paper are obtained from 
only the viewpoint of topology, without using any par- 
ticular models or hypotheses. These results are valid for 
all systems which support the existence of scroll wave. 

In this paper, we give a rigorous and general topologi- 
cal investigation of scroll wave. This work can be applied 
in practice. Here we give some discussions about how to 
connect our work to the phenomenology and mathemat- 
ical analysis of scroll wave. (i). The regular condition 
pni can be regarded as stability condition of scroll wave. 
When the solutions of Eq.([T|) is determined, the exact 
expression of pop can be work out. In this case, we can 
directly calculate the stability condition of scroll wave. 
This provide a direct approach to investigate the stabil- 
ity of scroll wave and how to control it. (ii). When the 
regular condition pTI)) fails, the scroll wave is unstable. 
In this case, the branch processes occur. The branch 
condition D{(f>/x) = predicts where and how these pro- 
cesses will occur, if the scroll wave solution of Eq.([T|) is 
work out. In laboratory, these branch processes can be 
produced by using the branch condition. This provide 
an experimental approach to test our work. (Hi) In this 
paper, we have predicted that the velocity of the vortex 
filaments is infinite when they are being annihilated or 
generated. The similar pheriomena has been obtained in 
the phase-ordering system[33. We expect that this phe- 
nomena will be observed for scroll wave in laboratory. 
(iv). For knotted vortex filaments, the Hopf invariant 



will play an important role in control the behaviours of 
scroll wave. During the branch processes, the Hopf in- 
variant will protect by topology. We also hope that this 
will be proved by experiment, (v). At last, we point out 
that Eq.((T]) determines the solution of Z, and of course it 
determines the regular condition and branch condition. 
The investigation presents in this paper is based on an 
important precondition that the scroll wave solutions can 
be worked out from Eq.((T]). 
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